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We have investigated the superconductivity discovered in boron (B)-doped diamonds by means 
of 11 B-NMR on heteroepitaxially grown (111) and (100) films. 11 B-NMR spectra for all of the films 
are identified to arise from the substitutional B(l) site as single occupation and lower symmetric 
B(2) site substituted as boron+hydrogen(B+H) complex, respectively. A clear evidence is presented 
that the effective carriers introduced by B(l) substitution are responsible for the superconductivity, 
whereas the charge neutral B(2) sites does not offer the carriers effectively. The result is also 
corroborated by the density of states deduced by 1/T\T measurement, indicating that the evolution 
of superconductivity is driven by the effective carrier introduced by substitution at B(l) site. 



In 2004, Ekimov and coworkers discovered the super- 
conduct ivity(SC) in heavily boron (B)-doped diamond 
synthesized by high-pressure and high-temperature tech- 
nique in very low carrier concentration (4-5 x 10 21 cm -3 ) 
[l| . Many theoretical studies have stressed the similarity 
to MgB2 with a high-T c value of 40 K where strong cou- 
pling of the holes at the top of the valence band with op- 
tical phonons plays an important role @, 0, Q • Recently 
angle-resolved-photoemission spectroscopy (ARPES) re- 
vealed that the holes on the diamond bands play an es- 
sential role in determining the metallic nature of heav- 
ily B-doped diamond superconductors [5j. Soon after the 
first discovery, it has been also reported on the B-doped 
diamond film synthesized by microwave plasma-assisted 
chemical vapor deposition (MPCVD) method [6[. One of 
the advantages of this method is that B concentration 
can artificially tune over a wide range, which enables 
us to elucidate how a superconducting transition tem- 
perature (T c ) depends on a doping level of boron (tib) 
in the homogeneously B-doped diamonds @, Partic- 
ularly, Umezawa et al. have performed the systematic 
study on homoepitaxially grown (111) and (100) films, 
and revealed that T c for the (111) films is by more than 
two times higher than for the (100) films despite of an 
equivalent B concentration for both films (see inset of 
FigHj)[8]. This remarkable difference might contain some 
hints to understand the mechanism of superconductivity 
and gives a further insight for searching higher-T c su- 
perconductivity in the B-doped diamond and the other 
related materials. 

In this letter, we report on a microscopic evidence 
for evolution of superconductivity in B-doped (111) and 
(100) diamond films, based on the analysis of n B-NMR 
spectra and n (l/Ti) for two predominant boron sites. 

NMR measurements were performed in five different 
diamond films; thick and thin (111) heteroepit axial films, 
and thick and thin (100) heteroepitaxial films, and thick 



polycrystalline film, deposited on appropriate substrates 
for each deposition condition by MPCVD method^, @]. 
Here the thickness of "thick" and "thin" films in this work 
are 50~100/im and ~3/im, respectively. The nominal B 
concentration nB,siMS is derived from secondary ion-mass 
spectroscopy (SIMS) measurement. T c of each sample 
is determined by an onset of diamagnetism under zero- 
field cooling, which coincides with an offset of zero resis- 
tance. NMR spectra were obtained by fast Fourier trans- 
form (FFT) technique at a constant magnetic field using 
a conventional phase-coherent-type NMR spectrometer. 
The Knight shift of n B-NMR was determined by the 
relative shift from the resonance frequency of B(OH)3 
(K[B(OH)s]« 20 ppm). The nuclear spin-lattice relax- 
ation rate 1/Xi were measured at a constant magnetic 
field of 2.5 and 10 T. 

Figures [U(a) and QJb) show n B-NMR spectra at 4.2 
K for the (lll)-thick and (100)-thin films, which are 
noteworthy that T c = 5 K for the (lll)-thick film is 
higher than T c = 3.4 K for the (100)-thin film, although 
^b,sims = 2.8 x 10 21 (cm -3 ) for the former is smaller 
than nB,siMS = 8-7 x 10 21 (cm -3 ) for the latter. In both 
films, two spectra overlap with a narrow peak around 
Af = with a linewidth of 5 kHz and broad ones in 
the range Af = —40 and +20 kHz, respectively. For 
the narrow peak, the Knight shift is approximately -10 
ppm, and for the broad one their shifts are distributed 
over a wider range, both of which do not depend on tem- 
perature in the range of 1.4 - 200 K. The narrow spec- 
trum originates from the most symmetric B site, i.e., the 
most likely substitutional position for the carbon site, the 
linewidth of which should be narrow due to the uniax- 
ial symmetry along each [111] direction at the carbon 
site in the diamond structure. We denote it as B(l) 
site hereafter. On the other hand, the broader spectrum 
ought to arise from the borons in lower local symme- 
try, such as boron+hydrogen (B+H) complex site, the 
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FIG. 2: (color online) (a) The broad peak becomes narrower 
in the higher field in all of the films. (b)The field depen- 
dence of the spectral edges is well reproduced by the case of 
vq ~1.8MHz as indicated by solid lines. It provides a clear 
evidence for presence of EFG only in B(2) site. 



FIG. 1: (color online) 11 B-NMR spectra for (a) (lll)-thick 
and (b) (100) -thin films, which are typical that T c for the for- 
mer is higher than that for the latter, although ns,siMS for 
the former is about three times smaller than that for the lat- 
ter. The calculated spectrum (red curve) for the (lll)-thick 
film is well reproduced by superposition of two components 
originated from B(l) and B(2) sites. 



interstitial B sites, B+B paired occupation sites, and so 
on. We identify this line-broadening as arising from ad- 
ditional nuclear quadrupole interaction derived from an 
electric field gradient (EFG) at n B nuclear site. It is 
discriminated as B(2) site. Generally, in case of n B nu- 
cleus (1=3/2) where the nuclear quadrupole interaction 
is sufficiently smaller than Zeeman interaction, the spec- 
trum for resonance peak(I=+l/2<^>-l/2) are broaden by 
the second order perturbation of EFG, which is given 
by A/(0) = Zv 2 Q sin 2 (9(1 - 9 cos 2 0)/(16j N H), where v Q 
is nuclear quadrupole frequency that is proportional to 
EFG, 7at is a nuclear gyromagnetic ratio and 6 is an an- 
gle between the principal axis of EFG and external field 
that should be taken as a random distribution in this 
material. The right and left edges given when = 90° 
and 6 = 41.8° follow the relation of +3i/%/(l&y N H) 
and —VqK^nH) against respectively. Actually, the 
spectral edges of the broad peak indicated by arrows in 
Fig. [21(a) become narrower in the higher field, the field 
dependence of which is well reproduced by the case of 
vq ^1.8MHz for (lll)-thin film, as indicated by solid 
lines in Fig. [21 b). This behavior has been seen in all of 
the films in our study, evidencing the presence of EFG 
only for B(2) site. Among those possible lower symmetric 
B sites, we consider that the B substitution for carbon 



site as B+H complex is most likely than others taking 
account of the MP-C VD p rocess using the mixed gas of 
CH4, (CH3VB and H2[l2[, as discussed also by ab initio 
calculation [uj [Til] - In this case, the presence of posi- 
tively charged hydrogen makes the spectrum broader in 
association with EFG as discussed above, and does not 
offer the carriers effectively since the positively charged 
hydrogen compensates the carrier doped by B substitu- 
tion. In fact, the recent SIMS measurement in these films 
has revealed the presence of higher hydrogen density for 
lower T c films [15]. In this context, the narrow and broad 
NMR spectra are identified as arising from the substitu- 
tional B site as single occupation (B(l) site) and as B+H 
complex (B(2) site), respectively, and are totally repro- 
duced by the composition of two different spectra from 
B(l) site and B(2) site with reasonable parameters [9j, as 
shown in Fig. [2(a). 

It is remarkable that the intensity of the B(l) site is 
predominantly observed for the (lll)-thick film, whereas 
not for the (100)-thin film. Since the NMR intensity is 
proportional to the number of B nuclei, we can estimate 
a fraction of B(l) site (/b(i)) against nominal B density, 
by means of evaluating a fraction of the narrow spectrum 
against an intensity integrated over a whole spectrum. 
The estimated /b(i)' s are about 49% for (lll)-thick film 
and 12% for (100)-thin film. With respect to the (100)- 
thick film (n B ,siMS = 3.8 x 10 21 cm- 3 , T c -5 K) and the 
polycrystalline films (nB,siMS = 4.4 x 10 21 cm -3 , T c ~5.2 
K), the NMR spectral shapes are somewhat similar to 
that of (lll)-thick film, as shown in Fig. [3j The surface 
observed by the scanning electron microscopy (SEM) are 
found to be not flat any longer due to the presence of 
the preferential {lll}-growth facets. The /b(i)' s are es ~ 
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timated to be ~33% and ~44% for the polycrystalline 
and the (lOO)-thick films, respectively. 
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FIG. 3: (color online) NMR spectra and the SEM images of 
surface for (100)-thick and polycrystalline films. These spec- 
tral shapes seem similar to that of (lll)-thick film, because 
both surfaces are dominated by the {111} growth facet. 

The B-concentration dependence of T c are quite impor- 
tant to unravel an origin for the onset of superconductiv- 
ity in the B-doped diamond. We propose that T c should 
be plotted as a function of the boron concentration at 
B(l) site (n B (i)) that is determined by nB,siMS x /b(i)? 
instead of a nominal concentration nB,siMS- As shown 
in Fig. 21 it has been found that T c 's in several dif- 
ferent grown films are plotted on a rather simple curve 
as a function of riB(i)- This scaling behavior presents 
clear evidence that effective carriers introduced by the 
B substitution as B(l) site are responsible for the on- 
set of superconductivity. On the other hand, the result 
also demonstrates that the B(2) site with B+H complex 
does not contribute to superconductivity because of the 
substitution as a neutral charged state. It is also cor- 
roborated that the effective carrier density estimated by 
ARPES[5],[H is also plotted on the extrapolated curve 
in Fig. 0J 

This result is also corroborated by measurements of 
the nuclear spin-lattice relaxation rate \jT\ for (111)- 
thick and (100)-thin films. As shown in Fig. Efa), l/T{T 
stays constant against T, obeying so-called "Korringa's 
relation" over a whole T-range, demonstrating that SC 
is realized under the Fermi liquid state for both films 
despite the vicinity of metal-insulator-transition(MIT) 
n c « 3 x 10 20 cm- 3 [14]. The relation of \jT{T oc N(E F ) 2 
enables us to deduce the density of states at the Fermi 
surface N(E^). As shown in Fig. E^b), it was found 
that (lll)-thick film with higher effective carrier den- 
sity has larger N(Ep) than (100)-thin film with lower 
carrier density does, suggesting that the enhancement of 
T c in (lll)-thick film primarily relates to the increase 
of N(Ep) that is induced by the effective carrier dop- 
ing by the substitution at B(l) site. Note that their T c s 
were not precisely predicted only by increase of N(Ep) 



FIG. 4: (color online)T c is plotted as a function of n B (i), in- 
stead of a nominal concentration tib,sims. A scaling behavior 
is found between T c and n B (i) corresponding to the effective 
carrier density. The crosses denote the carrier density esti- 
mated from ARPES 0, . 



in McMillan's formula [16j when assuming the same De- 
bye frequency and electron-phonon interaction, suggest- 
ing that it is necessary to consider the doping dependence 
of these parameters to describe its T c . 

Finally we comment on the result in the (lll)-thin 
film (n B ,siMS = 8.4 x 10 21 cm- 3 ) with highest-T c of 7 K 
in this work. Although T c is two times higher than that 
for (100)-thin film in spite of very similar nB,siMS? the 
apparent differences has not been observed in their spec- 
tral shape, as shown in Fig. EJa). It must be noted that 
the (lll)-thin films tend to exhibit broader SC transi- 
tion than the (100) films does(see Refs.0, 0, Hj]). Indeed, 
the diamagnetic response for this (lll)-thin film appears 
over the broad T-range between 2 K and 7 K as indi- 
cated in Figs. Efb) and [EJc). Recently a possible distri- 
bution of T c in the (lll)-thin films has been thoroughly 
investigated by grinding the (lll)-thin film from surface 
and/or bottom, and as a result, it is revealed that T c be- 
comes higher close to the surface by means of transport 
measurement [TBj . It suggests the spatially distribution of 
T c for the (lll)-thin films over a range in T c =2-7 K. The 
intensity of NMR spectra reflects the spatially- averaged 
B(l) and B(2) densities for all domains with different 
T c 's. In this situation, we tentatively supposed T c ~4 K 
as a spatially- averaged value, determined by a center of 
the broad peak in dx/dT vs T curve in Fig. Etc). Using 
the riB(i) ~ 1 x 10 21 cm -3 from the spectrum analysis, it 
is reasonably plotted on the same curve in Fig. [H sug- 
gesting the spatially- averaged T c in the (lll)-thin film 
seems to be related to n B (i) as well. In addition, \jT\T 
in (lll)-thin film is also similar to that in (100)-thin film, 
as shown in Fig. El suggesting that both N(E-p) and rism 
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FIG. 5: (color online) (a)T-dependence of l/TYT for (111) 
and (100) films. (b)N(E F ) deduced by y/(l/TiT) are plotted 
against effective carrier density n B (i)- (lll)-thick film with 
higher tib(i) has larger N(Ef) than (lOO)-thin film with lower 
tib (i) does, suggesting that the enhancement of T c in (111)- 
thick film relates to the increase of N(Ef). 



the hydrogen forming the B+H complex, but it has not 
been succeeded yet. 

In conclusion, the present 11 B-NMR study on the B- 
doped diamond revealed that borons are doped into the 
substitutional B(l) site as single occupation and as lower 
symmetric B(2) sites forming boron+hydrogen (B+H) 
complex. This result has enabled us to extract the scaling 
behavior between T c and the B(l) concentration 713(1)- 
A clear evidence is presented that the charge carriers in- 
troduced by B(l) substitution are responsible for the su- 
perconductivity, whereas the charge neutral B(2) sites 
substituted as B+H complex does not offer the carriers 
effectively. A route to make T c higher in the boron doped 
diamond is to dissociate the hydrogen at B(2) site by an- 
nealing, and to find how it can be possible to make the 
boron concentration at the substitutional B(l) site in- 
crease, keeping the diamond structure. This may also 
help developing diamond-based devices that make use of 
the unique properties of diamond. 
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are also similar with each other. 



(a) f=34.8B7MHz 
H=2.554T 

(11 1)-thin film T=4.2K 
n n = 8.4 x10 21 cm 3 



B(1) 
B(2) Ji 





(b) / /^-"~""*~* 

* 1 L^r ( 111 >- thin 

(100)-thin / 

/ (111)-thick 


(c) /?0 

(100)-thin 


\(111)-thlck 
I — \~^(111)-thin 







-80 -60 -40 -20 20 40 60 80 1 2 3 4 5 

A frequency ( kHz ) 

T(K) 



FIG. 6: (color online) (a) NMR spectra, (b) susceptibility 
and (c) its T- differential for (lll)-thin film, together with 
the others. The diamagnetism in (lll)-thin film gradually 
occurs over the broad T-range between 2 K and 7 K. 

However, there still remains a question why thin films 
includes more B(2) than thick film does. With respect 
to B+H complex, there are possible H positions such as 
bond-center, back-bond positions of boron-carbon, C2 V 
symmetry site [Tol. El|, B+H n complexes jl7j as suggested 
by the band calculations. Recently Oguchi calculated 
the EFG for these states and obtained a similar value 
of vq to our experiment for the case of the bond-center 
position of B-CpjJ. Either the relaxation process of the 
lattice constant in the thick films or a non-equilibrium 
process of MP-CVD method may control the formation 
of B+H complex. The annealing has been tried to expel 
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